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A transferable tight-binding potential has been constructed for heteroatomic systems containing carbon and hydrogen. The 
electronic degree of freedom is treated explicitly in this potential using a small set of transferable parameters which has been 
fitted to small hydrocarbons and radicals. Transferability to other higher hydrocarbons were tested by comparison with ab 
initio calculations and experimental data. The potential can correctly reproduce changes in the electronic configuration as a 
function of the local bonding geometry around each carbon atom. This type of potential is well suited for computer simulations 
of covalently bonded systems in both gas-phase and condensed-phase systems. 
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I. INTRODUCTION 

Computer simulation has become an important the- 
oretical tool in our understanding of complex molecu- 
lar systems. The two simulation techniques in statisti- 
cal mechanics - molecular dynamics and Monte Carlo 
- commonly used to explore the microscopic origin of 
macroscopic behaviors of many-body systems require as 
their input realistic microscopic interactions between the 
constituents. In this paper, we describe a general ap- 
proach for constructing transferable potentials for het- 
eroatomic systems based on tight-binding Hamiltonians. 
A potential for systems containing carbon and hydro- 
gen is reported here. It is part of an ongoing project 
to construct potentials for reactive systems of interest to 
organic chemistry, biochemistry and materials research. 

Systems containing carbon and hydrogen atoms hav- 
ing strong covalent bonds are interesting because of their 
central role in organic chemistry. Carbon atoms can as- 
sume three different hybridizations - sp, sp^ and sp^ - 
depending on the local and global chemical environment, 
and because of this carbon forms single, double, triple 
and conjugated bonds in compounds of many different 
geometries with tetrahedral, planar, linear, ring and cage 
structures. To understand these systems using computer 
simulations, it is necessary to first develop an interac- 
tion potential for carbon-hydrogen systems which could 
account for the complexities of the carbon-carbon and 
carbon-hydrogen covalent bonds. 

Pairwise potentials cannot describe covalently bonded 
hydrocarbons. A recent effort to construct potentials 
that can describe reactions in covalently bonded systems 
was first made by Tersoff who used an empirical bond- 
order type potential to describe primarily silicon-silicon 
and other covalent bonds. This method was later ex- 
tended to hydrocarbons by Brenner |^ who used a set 
of complex functional forms to parametrize the variation 
of the bonding geometry on each carbon atom based on 
its first and second nearest-neighbor carbon and hydro- 
gen atoms. The resulting potential involves a large set of 
parameters. 



In another recent approach, a potential based on the 
tight-binding Hamiltonian has been used to describe co- 
valent bonding in homoatomic systems containing either 
carbon or silicon ^ as well as adsorbed hydrogen 
on silicon surfaces po|-p^. This approach treats the 
electronic degrees of freedom explicitly using a small 
set of transferable parameters. Chemical bonding re- 
sults from filling the single-particle orbitals obtained 
from the diagonalization of the tight-binding Hamilto- 
nian. This alternative approach to constructing inter- 
atomic potentials has many advantages over conventional 
parametrization schemes. First, because of its quantum 
mechanical origin, this potential has the requisite func- 
tional form that is able to describe bond breakage and 
bond formation. Second, only a small set of parameters 
is necessary to parametrize a large class of molecular in- 
teractions. Third, the potential is highly transferable, 
i.e. the parametrization is obtained by fitting to only a 
very small set of chemical species (typically four or five), 
but the potential is extendible to very large systems that 
may be substantially more complex. 

In this letter, we report a new transferable minimum- 
basis tight-binding (TB) potential for hydrocarbons. We 
believe this is one of the first successful applications of the 
tight-binding potential to describe covalent bonding and 
chemical reactions in heteroatomic systems. Our results 
on hydrocarbons indicate that using this approach we are 
able to generate potential energy surfaces which compare 
well with high-level ab initio calculations. This potential 
is well suited for dynamical simulations because of its 
simplicity and its ability to account for a large variety of 
bonding configurations. 



II. METHOD 

Hydrocarbons have strong covalent two-centered direc- 
tional bonds. To describe the hybridization and the di- 
rectional bonding character in hydrocarbons, we adopt 
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a semi-empirical minimum basis tight-binding Hamilto- 
nian. The tight-binding Hamiltonian, which considers 
only the 2s and 2p valence electrons of carbon and the 
Is electron in hydrogen, is given by 



H' 



TB 



(1) 



where i and j label the atoms and a and (3 label the 
atomic orbitals. is the atomic orbital energy of atom i 
and orbital a. ^ is the hopping matrix element be- 
tween atomic orbital a on atom i and (3 on atom j. 
Eq. (P is the conventional tight-binding Hamiltonian of- 
ten used for sohds 

In our system, there are four types of hopping matrix 
elements between carbon atoms: t^S, tj^^, t^^, i^^, 
and two types of hopping matrix elements between car- 
bon and hydrogen: i^j^, t^p„. The detailed forms of the 
hopping matrix elements will be discussed later. Note 
that since there is no electron-electron interactions in 
Eq. (Q), 7?TB can be exactly diagonalized as a noninter- 
acting system, yielding the molecular orbitals. The elec- 
tronic configuration is then defined by putting the req- 
uisite number of electrons into the one-electron orbitals 
subject to the exclusion principle. 

The cohesive energy of the system is defined as 



Ecoh — E, 



Val 



E 



Ecorci^ij) 
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where -Evai gives the total energy due to the valence elec- 
trons, the core repulsion energy -E^^re comes from the 
screened ion-ion interactions between atoms i and j, and 
the atomic energy E^^^^ is the reference energy of the 
isolated atom i in the dissociation limit. 

The valence energy .Evai comes from two terms: 



(3) 



The first contribution to i?vai comes from the molecular 
orbital energies, summed over all orbitals resulting from 
diagonalizing the tight-binding Hamiltonian. n-y is the 
occupation number of orbital 7, and is the single par- 
ticle energy of orbital 7. In other words, the first term 
in Eq. (|^) is simply the contribution of H^b to the total 
energy. The second term in Eq. (^) places a restriction 
on the electrons occupying the same molecular orbitals. 
This [/-term describes empirically the extra Coulomb re- 
pulsion and correlation energy that arise when two elec- 
trons are in the same molecular orbital 7. This term 
contributes only when the molecular orbital is fully oc- 
cupied, i. e., the occupational number is 2. It is similar 
in spirit to the U-tevm in the Hubbard model ||l^ and is 
included in an attempt to describe electron correlation in 
an approximate way. Because of the /7-term {U > 0), an 
energy penalty is paid each time a pair of electrons (of 
different spins) go into the same orbital. Therefore it is 
possible in some circumstances that total energy is made 



more favorable by exciting an electron from a doubly oc- 
cupied level to an unoccupied orbital higher in energy 
to avoid paying the penalty. Consequently, the valence 
energy ii'vai has to be minimized with respect to the oc- 
cupation numbers {ri'y}. 

For the tight-binding Hamiltonian, we use the Slater- 
Koster parameterization scheme for the electronic 
hopping matrix elements. We adopt scaling forms |^ for 
the distance dependence of the hopping matrix elements 
as well as for the core repulsive interactions. 



ia.,l3{r) = ta,l3 
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(4) 
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where rj and are cutoff distances for the hopping ma- 
trix elements and repulsive interactions, and the coeffi- 
cients and power exponents ria, nc, and m^, m^, rric 
determine the general shape and sharpness of the cut- 
off functions. For the same pair of atoms, the reference 
bond distance rg is the same for both the hopping ma- 
trix element and the core repulsion. Unlike the original 
formulation ||^, however, we do not demand the hopping 
parameters to have the same power exponent as the 
core repulsion rric- 

The tight-binding Hamiltonian explicitly considers the 
contribution of the electrons to the total energy of the 
system. It is a natural extension of the often used 
extended Hiickel model in organic chemistry [ p^ . In 
strongly covalently bonded systems, the molecular or- 
bitals are localized and the chemical bonds are direc- 
tional. Consequently, the overlap of the wavefunctions 
are mainly along the line connecting two atoms, and 
two-center integrals dominate the electron-electron inter- 
actions, while three-center and four-center integrals are 
negligible. By parameterization of the hopping matrix 
elements, we have taken into account the two-center in- 
teractions. Correlation effects, which are very important 
in certain circumstances, are completely missing in the 
conventional molecular orbital theory. Some correlation 
effects have been reintroduced into our potential via the 
Hubbard-like C/-term in the Hamiltonian. Although the 
?7-term seems to have been implemented in an ad hoc 
manner, it plays a central role in determine the occupa- 
tion of the molecular orbitals. In particular, it guaran- 
tees that molecules dissociate correctly without anoma- 
lous charges in the separated atom limit. For example, a 
CH radical in the tight-binding Hamiltonian without [/ 
will dissociate to give the unphysical products H^ + C+ . 
In addition to producing the total energy of the system, 
the tight-binding Hamiltonian also generates the single 
particle energy states. These single-electron energies are 
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useful for an approximate picture of the electronic den- 
sity of states. 

There are two parts to the potential: the carbon- 
hydrogen interactions and the carbon-carbon interac- 
tions. We will briefly describe how the parameters of 
these two parts of the potential were obtained. 



A. Carbon-Hydrogen Interactions 

For the carbon-hydrogen interactions, the hopping 
matrix elements and repulsion parameters at the refer- 
ence bond distance tq, tf^{ro), tf^^ri'^o), and Ef^^{ro) 
were first fitted to the experimental energies for CH, 
CH2, CH3 and CH4, taken from Ref. [Q. Since the 
C-H bondlengths in these four CH species are similar, 
we set the reference bond distance to vq — 1.09A. Care 
was taken to ensure that the resulting hopping matrix 
elements would reproduce a planar equilibrium structure 
for CH3. Then the power exponents Ha for the hopping 
matrix elements, the power exponent rua for the core re- 
pulsion and the cutoff functions were chosen to give the 
correct experimental equilibrium bond distance of CH4 
as well as the correct breathing mode frequency of CH4. 
The resulting parameters are listed in Table ^. 



B. Carbon-Carbon Interactions 

There are currently available several sets of carbon- 
carbon tight-binding parameters 17 1. For this work, 
we have developed a new set. For the carbon-carbon 
interactions, the four hopping matrix elements ^^^(ro), 
t^peri'i^o), tppA'T'o), and i^p^(ro), and the repulsion param- 
eter i?TO?c(^o) were fitted to the carbon-carbon bond en- 
ergies of C2, C2H2, C2H4, and C2H6 at ro = 1.312A. 
The bond energies of C2 and C2H4 were taken from ex- 
periments jl^jl^, whereas the bond energies of C2H2 
and C2H6 at ro = 1.312A were deduced from the re- 
spective bond dissociation energies at their equilibrium 
distances and the corresponding force constants of the 
carbon-carbon bonds. Then the power exponents for 
the hopping matrix elements, the power exponents for 
the core repulsion and the cutoff functions were fitted 
to the equilibrium carbon-carbon bondlengths and force 
constants in C2, C2H2, C2H4, and C2H6. The result- 
ing carbon-carbon interaction parameters are listed in 
Table |. 



TABLE I. Parameters for the tight-binding potential. 



III. RESULTS 

In this section, we will analyze specific features of our 
tight-binding (TB) hydrocarbon potential. A compar- 
ison of the TB potential to ab initio calculations and 
experimental data is summarized in Tables || to VI , 
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A. Structures 

In Table |l^, the equilibrium geometries of small hy- 
drocarbons from the TB potential are compared to ah 
initio calculations and experimental results [ p^ . 

Excellent agreement with the experimental values was 
found. For the species with only carbon-hydrogen inter- 
actions, i.e. CH4 and CH3, our semiempirical potential 
reproduces the correct equilibrium geometries and the 
correct carbon-hydrogen bondlengths to within 0.01 A. 
For the molecules involving not only carbon-hydrogen in- 
teraction but also carbon-carbon interactions, i.e. C2H2, 
C2H4, C2H6, CeHg, excellent agreement with experimen- 
tal results is also achieved. The carbon-hydrogen bonds 
in these four hydrocarbons are only slightly longer com- 
pared to experimental data and the experimental trend 
that the carbon-hydrogen bondlength increases from 
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TABLE II. Tight-binding, ab initio and experimental equilibrium geometries for simple hydrocarbons. All ab initio results 
are derived from calculations using a 6-31G** basis set unless otherwise indicated. All bondlengthes are in A. 



Molecule 


Point group 


Feature 


TB 


HF'' 


MP4 


Expt.'' 


CH3 


Dsn 


r(C-H) 


1.079 


1.075 


1.078 


1.079 


CH4 


Td 


r(C-H) 


1.094 


1.084 


1.094 


1.094 


C2H2 


Doah 


r(C-H) 


1.066 


1.057 


1.06 


1.058 


(triple) 




r(C-C) 


1.183 


1.186 


1.22 


1.208 


C2H4 


D2h 


r(C-H) 


1.094 


1.076 


1.08 


1.0836 


(double) 




r(C-C) 


1.341 


1.316 


1.34 


1.337 






ZCCH 


122.8° 


121.75° 


121.5° 


121.25° 


C2H6 




r(C-H) 


1.104 


1.086 


1.09 


1.091 


(single) 




r(C-C) 


1.546 


1.527 


1.53 


1.536 






ZCCH 


110.8° 


111.2° 


111.0° 


110.90° 






r(C-H) 


1.095 


1.08° 


1.09° 


1.084 


(conjugated) 




r(C-C) 


1.428 


1.38° 


1.39° 


1.397 



^ Ref. jigf 
° 3-21G. 



and |l|. 



sp to sp'^ hybridization is nicely reproduced. Our TB 
potential also reproduces correctly the trend that the 
carbon-carbon bondlength decreases from C2H6 to C2H6 
in the order single > conjugated > double > triple bonds. 
The carbon-carbon bondlengths in single and double 
bonds are within O.OlA from the experimental values, 
while conjugated and triple CC bonds are within O.OSA 
from the experimental data. The bond angles are also in 
excellent agreement with experimental measurements. 



B. Vibrational Frequencies 

Next we examine the properties of the potential en- 
ergy surface away from the minimum 



In Table g the 
vibrational frequencies are tabulated and compared to 
experimental and ab initio results ||lj,|l^ . Note that the 
fitting of the parameters involves only the C-H and C-C 
stretching mode frequencies, while all other frequencies 
are derived from the TB potential, and they reproduced 
the vibrational spectra remarkably well. Strong anhar- 
monic effects are observed in small hydrocarbons, and 
thes have been well established theoretically and ex- 
perimentally plf . 

The vibrational frequencies for CH3 and CH4 are all 
within 5% of the MP2 calculations except for the defor- 
mation mode in CH4 which is within 10%. The slightly 
larger apparent disagreement with the experimental data 
on the umbrella mode of the CH3 is due to strong anhar- 
monic effects. Figure Q shows the TB potential energy as 
a function of the C-H bondlength for CH4 constrained 
in the tetrahedral geometry. Clearly, the TB predictions 
agree well with MP4 results even for large deviations from 
the equilibrium structure. 

In the case of C2H2, the CC stretching frequency is 7% 
larger than the MP2 result. Other high energy CH modes 
are in excellent agreement with the MP2 results. The 



two low-energy bending modes seem stiffer in our empir- 
ical potential than the MP2 results, but the discrepancy 
with experimental data can be largely attributed to an- 
harmonicity. In the case of C2H4, all modes agrees well 
with the MP2 results, even the low energy ones. The 
CC stretching frequency is correct to within 3% com- 
pared to the experimental value. In the case of C2H6, 
in general all modes agrees well with the MP2 results. 
The CC stretch mode is 10% higher than MP2. The 
most noticeable discrepancy is the zero frequency pre- 
dicted by the TB potential for the torsional motion be- 
cause there is no barrier for the torsional motion in the 
TB potential for C2H6. This lack of a torsional barrier 
is a direct consequence of the symmetry of C2H6 within 
the minimum-basis TB Hamiltonian. The energy barrier 
in the torsional motion arises from the steric interaction 
between the third nearest-neighbor hydrogens and this 
component is currently unaccounted for by the present 
model. 



C. Energetics 

In chemical reactions, bonds are broken and new ones 
are reformed. The reaction energy will depend crucially 
on the energetics of different species involved in the re- 
action. In Table IV and 0, TB predictions of the en- 
ergetics involving the breakage of carbon-hydrogen and 
carbon-carbon bonds under different situations are com- 
pared to ab initio calculations and experimental data 
Mpl. In Table 0, 



the energetics of the sequen- 
tial removal of hydrogen atoms from methane are sum- 
marized. When each hydrogen is being removed, the 
remaining carbon-hydrogen bonds are allowed to re- 
lax to the minimum energy configuration. During the 
bond breaking process, we start from a stable closed- 
shell species CH4. Removal of one hydrogen pro- 
duces a planar CH3 radical with one unpaired electron. 
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TABLE III. Vibrational frequencies (in cm~^) for 
methane (CH4), the methyl radical (CH3), acetylene (C2H2), 
ethylene (C2H4), and ethane (CaHg), 



Molecule 


Symmetry 


Description 


TB 


MP2'' 


Expt.'^ 


CH4 


Oi 


s-stretch 


3162 


3115 


2917 




e 


d-deform 


1690 


1649 


1534 




t2 


d-stretch 


3252 


3257 


3019 






d-deform 


1570 


1418 


1306 


CH3 


a-i 


s-stretch 


3207 


3208 


3004 




02 


umbrella mode 


411 


404 


606 




e 


d-stretch 


3419 


3398 


3160 






scissor 


1552 


1478 


1288 


C2H2 


-t 


CH stretch 


3546 


3593 


3374 






CC stretch 


2146 


2006 


1974 






CH stretch 


3355 


3516 


3289 






CH bend 


811 


444 


612 






CH bend 


897 


783 


730 


C2H4 


a\g 


CH2 s-stretch 


3272 


3231 


3026 






CC stretch 


1680 


1724 


1623 






CH2 scissor 


1561 


1425 


1342 






CH2 twist 


1158 


1083 


1023 




big 


CH2 a-stretch 


3305 


3297 


3103 






CH2 rock 


1398 


1265 


1236 




blu 


CH2 wag 


1102 


980 


949 




b2g 


CH2 wag 


1084 


931 


943 




b2u 


CH2 a-stretch 


3344 


3323 


3106 






CH2 rock 


899 


873 


826 




bsu 


CH2 s-stretch 


3221 


3222 


2989 






CH2 scissor 


1651 


1523 


1444 


C2H6 


aig 


CH3 s-stretch 


3113 


3086 


2954 






CH3 s-deform 


1633 


1493 


1388 






CC stretch 


1157 


1040 


995 






torsion 





452 


289 




a2u 


CH3 s-stretch 


3148 


3104 


2986 






CH3 s-deform 


1614 


1494 


1379 




eg 


CH3 d-stretch 


3168 


3228 


2969 






CH3 d-deform 


1639 


1520 


1468 






CH3 rock 


1336 


1264 


1190 






CH3 d-stretch 


3201 


3215 


2985 






CH3 d-deform 


1621 


1604 


1469 






CH3 rock 


911 


783 


822 



^ 6-31G* MP2 calculations from ^ 
''References [H and itll. 



TABLE IV. Calculated and experimental bond energies 
(in kcal) for sequential removals of hydrogen from methane. 



Theory 


CH4 


CH3 


CH2 


CH 


TB 


107 


119 


106 


85 


HE'' 


87 


88 


101 


55 


MP2'' 


109 


110 


109 


73 


MP4"- 


110 


112 


107 


76 


Expt.'' 


112 


116 


107 


84 




1.5 



6-31G** calculations from Refs. idl and 18 



Ref. 16 



R(C-H) (A) 



FIG. 1. The cohesive energy of methane as a function of 
the carbon-hydrogen bondlength (constrained to maintain its 
Td symmetry): from MP4 6-31G** (circles), from TB (solid 
line). The MP4 results are shifted to the experimental cohe- 
sive energy. 

Removal of another hydrogen atom from CH3 produces 
a bent CH2 vifith two unpaired electrons. Finally, the 
removal of a third hydrogen produces CH with three un- 
paired electrons. The sequence of products obtained un- 
der sequential removal of hydrogen from CH4 involves 
a number of very different hybridizations, and hence 
the resulting electronic configurations and geometries are 
rather diverse. Table [V shows that our semi-empirical 
TB potential can reproduce the energetics quite satis- 
factorily. The errors in the sequential carbon-hydrogen 
dissociation energies are within 5 kcal from the experi- 
mental data, and the results are much better than ab ini- 
tio Hartree-Fock calculations. The order of the carbon- 
hydrogen bond strengths is also correctly reproduced. 

The carbon-carbon bond dissociation energies in 
C2H6, C2H4 and C2H2 are shown in Table^. In the case 
of C2H6, a single carbon-carbon bond is broken and the 
system dissociates from a closed-shell structure into two 
open-shell methyl radicals. The double bond in C2H4 is 
harder to break than the single bond and the end prod- 
uct is two open shell CH2 species. The triple bond in 
C2H2 contains the most energy among the three. Our 
empirical potential reproduces the three different bond 
energies excellently. 

The diversity of organic compounds relies on the abil- 
ity of carbon atoms to assume different hybridizations in 
order to form single, double and triple bonds. During 
chemical reactions, carbon atoms will rehybridize from 
one configuration to another continuously. A success- 
ful potential should be able to describe such processes 
closely. It will therefore be a more stringent test for 
our TB potential to see if it will reproduce the cor- 
rect rehybridization behavior along the reaction path. 
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100 - 



1.0 1,5 

R(CH3-H) (X) 

FIG. 2. The calculated (a) cohesive energy and (b) HCH3 
angle change due to the variation in the hybridization by grad- 
ually pulling one hydrogen atom out from methane: from 
MP4 6-31G** (circles), from TB (sohd line). 

In Fig. |^(a), we show the potential energy curve for the 
removal of one hydrogen from the tetrahedral CH4 to 
yield a planar CH3. In the process of pulling one hydro- 
gen out of the tetrahedron, the carbon atom reconfigures 
itself from sp^ to sp^ accordingly to open up the umbrella 
formed by the other three hydrogen atoms until a planar 
geometry is attained. In Figure ||(b), the HCH3 angle is 
plotted as a function of the H-CH3 distance during the 
rehybridization of the carbon atom. The predictions of 
the TB potential agree well with ah initio MP4 calcula- 
tions. 



D. Transferability 

By far, the most significant feature of the TB potential 
is its high degree of transferability. The TB potential 
has been fitted only to small molecules and radicals. In 
this section, we will demonstrate the transferability of 
the TB potential by comparing the TB predictions to 
experimental results for very much larger system. 

Extended structures can be produced by carbon- 
carbon bonds. One important example is the cage 
structure of Cgo. There are 12 pentagons and 20 
hexagons in Ceo and all 60 atoms are equivalent. 



TABLE V. Calculated and experimental carbon-carbon 
bond energies (in kcal) for single, double and triple bonds. 

H3C — CH3 H2C = CH2 HC = CH 



TB 94 177 235 

HF^' 69 119 167 

MP2'' 100 177 246 

MP4'' 98 177 236 

Expt.'' 97 179 236 



6-31G**. 
^ Refs. (l|l and (l|. 

TABLE VL Comparison of tight-binding predictions for 
atomization energies (in eV) of straight-chain alkanes. 



n Experimental TB 

1 18.22 18.13 

2 30.90 31.03 

3 43.3'' 43.90 

4 56.2*' 56.78 

5 69.0*= 69.65 

6 81. 8*^ 82.52 



"Ref. ||. 

'^Estimates corrected for zero-point energies. 

Because of the spherical cage structure, carbon atoms do 
not have either pure sp'^ or sp^ bonding. The chemical 
bonds in the pentagons, which are often called "single 
bonds", are 0.04 A longer than the remaining "double 
bonds" in the hexagons. This difference between the dou- 
ble and single bond distances is very nicely reproduced by 
the TB potential, although the absolute value predicted 
by the TB potential for a single bond is 0.05 A too long. 
The TB potential also predicts an atomization energy of 
7.12 eV per carbon atom, which compares very favorably 
with the experimental result of 7.02 eV per carbon atom 
11,11 . 

Another example of an extended structure produced by 
carbon-carbon bonds is the sheet structure of graphite. 
The experimental atomization energy of graphite is 
7.37 eV per carbon atom [p^ . Ignoring interplanar in- 
teractions, our TB potential predicts 7.40 eV per carbon 
atom. 

The benzene molecule which has conjugated CC bonds 
also presents an interesting test for the TB potential. 
The experimental atomization energy is 59.56 eV while 
the TB potential predicts 59.72 eV. This yields a C-C 
bond energy of 145 kcal, which is between the CC single 
and double bond energies given in Table 

For the first six straight-chain alkanes C„Il2n+2, we 
have listed the TB predictions for the atomization energy 
m Table 0. The agreement with experimental values 
is almost quantitative. The TB potential predicts an 
average increment of 12.9 eV per additional CH2 unit, 
which is very close to the experimental value 12.8 eV. 

Finally, we examine one cyclic alkane. The lowest en- 
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ergy conformation of cyclohexane is the chair form. Ex- 
perimentaUy, the chair form is more stable than the boat 
form by 5.5 kcaL This difference between the chair and 
the boat forms is due to steric hindrance which as dis- 
cussed above is unaccounted for by the current model. 
Therefore, in the present form, the TB potential predicts 
essentially no energy difference between them. The bar- 
rier between these two states is predicted to be about 4.8 
kcal/mol (with the transition state assumed to be the 
totally relaxed planar structure). Experimentally, this 
barrier is about 10.8 kcal/mol. 

IV. CONCLUSION 

We have described a transferable semi-empirical min- 
imum basis tight-binding Hamiltonian for treating reac- 
tions involving hydrocarbons. This potential treats the 
electronic degrees of freedom explicitly, describing the 
interatomic interactions using a small set of transferable 
parameters. Chemical bonding results from the diagonal- 
ization of the tight-binding Hamiltonian and then the fill- 
ing of the one-particle energy levels with electrons. This 
approach can be systematically extended to include more 
elements across the periodic table. Our results on hydro- 
carbons show that the electronic configuration depends 
sensitively on the local bonding configurations around 
each atom. The rehybridization of a carbon atom from 
sp"^ to sp^ under a change of geometry from a tetrahedral 
CH4 to a planar CH3 radical can be followed closely and 
agrees well with high level ab initio calculations. The 
energetics and topologies of all carbon-carbon bonds — 
single, double, triple and conjugated bonds — can be 
reproduced accurately. Transferability to substantially 
larger systems have been demonstrated. Due to its sim- 
plicity and its ability to describe a large variety of bond- 
ing topologies, our approach will be useful for dynamical 
simulations of reactive systems in both the gas phase and 
the condensed phase. 
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